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ABSTRACT 

The logarithmic potential is of great interest and relev ance in the study of the 
dynam ics of galaxies. Some sm all corrections to th e work of IContopoulos fc Seimenisl 
(|1990D who used the method of lPrendergastl (|1982l ) to find periodic orbits and bifurca- 
tions within such a potential are presented. The solution of the orbital radial equation 
for the purely radial logarithmic pot ential is then considered using the p-cllipsc (pro- 
cessing ellipse) method pioneered by IStruckl (|2006f ) . This differential orbital equation 
is a special case of the generalized Burgers equation. The apsidal angle is also deter- 
mined, both numerically as well as analytically by means of the Lambert W and the 
Polylogarithm functions. The use of these functions in computing the gravitational 
lensing produced by logarithmic potentials is discussed. 
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1 INTRODUCTION 

The logarithmic potential has great interest in connection 
with the dyna mics of elliptical g alaxies and galactic halos. 
Introduced by iRichstong (|l98Cf ) to model stellar systems 
with concentric axisymmetric oblate spheroidal potential 
surfaces, it is one of the few axisymmetric galactic potentials 
with an equally simple mass distribution function. As a re- 
sult it has been studied extensive ly (e.g. iBinnev fc Spergell 
d 19821) :lBinnev fc Tremaine] l|l987l )). 

Richstonel (|l982t )~ did an extensive survey of orbits 



tials have been exam ined analytically JBelmonte et al.l 20071 : 
iPucacco et al.ll2008l ) and numerically I Magnenaj|l982 ). 



within scale-free logarithmic potentials (that is, with zero 
core radii). The effect of core radius and the presence 
or absence of a c entral mass have been exam i ned b y 
Gerhard fc Binnevl il 19851): iPfenniger fc de Zeeuwl (1 19891) : 



Miralda-Escude fc Schwarzschildl (|l989l ). lEvand (|l993h 



amined the axisymmetric case of galaxies embedded in ex- 
tended dark matter halos. iLees fc Schwarzschildl { 19921) ex - 
amined triaxial halo models. Karanis fc Caranicolasl ( 200ll ) 
examined how the core radius and the angular momentum 
are related to transitions from regular motion to chaos in log 
potentials. iTouma fc Tremaind (|l997l ) developed a symplec- 
tic map to study the dynamics of orbits in non-spherical 
potentials, with particular emphasis on the logarithmic 
potential. Periodic orbits in triaxial logarithmic poten- 



Email: valluri@uwo.ca 



Beyond galactic dynamics, the potential also has ap- 
plications to the problem of gravitational lensing. Beyond 
astrophysics, applications of the logarithmic potential oc- 
cur in the solutio n of planar boundary value problems in 
potential theory (|Evansl Il927l ) and with boundary value 
problems of analytic function theo ry. In elementary parti- 
cle physics, iQuigg fc Rosnerl \ 19771 ) use the logarithmic po- 
tential to show that the quarkonium level spacings are in- 
dependent of quark mass, in the non-relativistic limit. In 
this paper, the analysis of IContopoulos fc Seimenisl (|l990l ) 
( CS) is re-examined. CS applied the analytical techniques 
of iPrenderg ast (1983) to find approximate solutions to the 
equations of motion for particles moving within a logarith- 
mic potential. The Prendergast method was introduced to 
approximate some complex differential equations, such as 
the Duffing equation, and new applications for this method 
are still being found today. We then elaborate on the work 
of CS, turning our attention to the radial orbital equation, 
using the nonlinear Burgers equation to determine an ap- 
proximate analytic solution from which the apsidal angle is 
determined. It is also determined by finding the roots of the 
Lambert W and the polylogarithmic function. 

In sections[2]and|31 the Prendergast Method is revisited. 
We performed a thorough stud y of the pioneering work of 
IContopoulos fc Seimenisl (|l990l ) and present a slight elabora- 
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tion as well as a few mino r correc t ions t o their equations. In 
section[4] we briefly studv lstruckl |2006l Vs p-ellipse (precess- 
ing ellipse), introduced in his fine work on precessing orbits 
in a variety of power-law potentials, some shallower than the 
1/r Keplerian one. These potentials include the logarithmic 
potential of zero as well as nonzero core softening length. We 
present an integrable equation that provides us with values 
for the apsidal angles of the orbits considered. In section [5] 
we discuss the deflection of light in a logarithmic potential 
and gravitational lensing. Finally, section [6] summarizes our 
conclusions and any further work to be considered. 



2 UNPERTURBED SOLUTIONS 

Here we begin by re-establishing the results of CS with some 
minor corrections, before going on to use these solutions in 
subsequent sections. Where alterations to their values are 
given, they are indicated by asterisks. 

Following the notation of CS, our expression for the 
logarithmic potential is 



V(x,y) = In 



x 2 + ^ 

x +(72 



c- 



(1) 



where C is the core radius and U describes the elliptic- 
ity of the potential. Though of mathematical interest over 
a wider range of parameters, models with U > 1.08 or 
U < l/V% = 0.707 are unphysic al in that they require nega- 
tive mass densities l|Evanslll993h . As a result, only values of 
0.707 < U < 1.08 are of interest to galactic dynamics. The 
CS method begins by finding a solution for arbitrary values 
of U in the one-dimensional case (y = y = 0), and adding 
the motion in the second dimension as a perturbation. 

By finding the derivative of Eq.[T]and introducing it in 
the relevant second order orbital differential equation, it is 
possible to develop two equations of motion-one for the x 
component, the other for the y component: 



x + 



2U 2 x 



V + 



U 2 x 2 + y 2 + C 2 U 2 
ly 



U 2 x 2 +y 2 + C 2 U 2 



(2) 
(*) (3) 



where the (*) indicates the equation contains a correction 
to a typo in CS's original. 

The subsequent s olution is developed using the method 
of IPrendereastl [jl982h . Developed for second-order nonlin- 
ear ordinary differential equations, Prendergast applied the 
technique to the van der Pol oscillator and Duffing's equa- 
tion. CS applied it to the orbital equation in the logarithmic 
potential. 

The method begins by assuming a solution for x and y 
of the following form: 

N M ^ 

x = y = (4) 

D D w 

where N, M and D are Fourier series of the form 
N = N k cos(kujt), 

k— odd 

D — 1+ Acos(Zwi) (5) 



and which are truncated at the appropriate order. In the 
unperturbed one-dimensional case, y = y = and M = 0. 

In determining the solution for x, we introduce the ex- 
pansion of [4] into Eq. [2j and solve for a new equation of 
motion, 



(N"D 2 - 2N'D'D - ND"D + 2ND' 2 ) x 
(U 2 N 2 + M 2 + C 2 U 2 D 2 ) + 2U 2 ND 4 = 0. 



(6) 



A solution is essayed of the form 

TV = Acos(u)t); D = 1 + B cos(2ojt); (7) 

with constants A, B and lu to be determined, though we 
require B ^ for a non-trivial rational approximation. 

Finally, we introduce the proposed solutions[7]into Eq.[S] 
and set equal to zero the coefficients of cos (art) and cos(3oj£). 
This gives us two equations 

cj 2 fci + 0.75B 4 + 3B 3 + 6B 2 + 4B + 2 = 

u?k 2 + 0.5S 4 + 3B 3 + SB 2 + 4B = (*) (8) 

where fci and hi are given below; 

ki = 3.5625A 2 B 2 + 7B 2 C 2 + 2.125C 2 B 4 - (9) 

2BC 2 + 3.5C 2 B 3 - 0.75A 2 - C 2 , 
k 2 = 0.1875A 2 B 2 + 1.25A 2 B- 0.25A 2 + (10) 

6.5C 2 B 3 - 3.5B 2 C 2 + 2BC 2 + 0.25C* 2 B 4 . (*) 

The third equation needed to determine A, B and uj is given 
by the initial condition 



(1 + B)x Q - A = 



(11) 



where xo = x{t = 0). 

We now solve these equations for A, B and ui with 
the given values of xq. The solutions, as well as all mathe- 
matical manipulations presented in this paper unless other- 
wise mentioned, were determined using the software package 
Maple 15. The solutions have cj 2 > and are presented in 
Tabled] 

We note that for motion solely in the x-direction the 
value of U is irrelevant, and it appears neither in Eq. [8] nor 
in the initial conditions. 

An example of an unperturbed solution is displayed in 
Figure [T] with values corresponding to line 17 in Table [T] 



3 PERTURBED SOLUTIONS 

Purely radial orbits such as those of Section[2]are unlikely in 
practice. Here the search is for solutions to the motion where 
the y-component of the motion is close to the unperturbed 
motion discussed previously. In this case, M is no longer 
identically zero and we look for solutions of the form 



M — Mo + 5M; D = D + SD; 



N = N + SN (12) 



where the subscript indicates the unperturbed solution. 
The next st ep is to solve the differential equ ation, introduced 
as Eq. 8 in IContopoulos fc Seimerua (|l990h 

(U 2 N 2 D 2 + C 2 U 2 D A )5M" - 

2 7T 2 n 3 rV\ 



(2U 2 N 2 DD' + 2C 2 U 2 D S D')6M' + 



(2U 2 N 2 D' 2 - U 2 N 2 DD" 



C 2 U 2 D 3 D" + 



2C 2 U 2 D 2 D' 2 + 2D 4 )5M = 0. 



(13) 
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xo A B u 



1 


0.001 


0.001000006 


6.2497E-6 


14.142 


2 


0.01 


0.01000621694 


0.0006216935078 


14.089 


3 


0.02 


0.02004895506 


0.002447752803 


13.935 


i 


0.03 


0.03016097774 


0.005365924678 


13.689 


5 


0.04 


0.04036817727 


0.00920443171 


13.368 


6 


0.05 


0.05068761968 


0.01375239358 


12.989 


7 


0.06 


0.06112700808 


0.01878346808 


12.569 


8 


0.07 


0.07168549791 


0.02407854158* 


12.125 






(0.024074) 






9 


0.08 


0.08235549842 


0.02944373* 


11.668 






(0.029436) 






10 


0.09 


0.09312493232* 


0.03472147* 


11.211 






(0.093124) 


(0.034708) 




11 


0.1 


0.1039794453* 


0.039794453* 


10.761 






(0.103977) 


(0.039775) 




12 


0.125 


0.13139283* 


0.05114264* 


9.696 






(0.131388) 


(0.051102) 




13 


0.15 


0.1590498684* 


0.060332456* 


8.747 






(0.159040) 


(0.060266) 




11 


0.175 


0.1868185967* 


0.067534838* 


7.921 






(0.186803) 


(0.067443) 




15 


0.2 


0.2146246114* 


0.073123057* 


7.210* 






(0.214602) 


(0.073008) 


(7.209) 


16 


0.225 


0.2424301778* 


0.077467457* 


6.597 






(0.242400) 


(0.077332) 




17 


0.25 


0.27021796* 


0.08087184* 


6.068 






(0.270180) 


(0.080719) 





Table 1. A comparison of our results with those of Contopolous 
and Seimenis (1990). Where a corrected value appears, the origi- 
nal value appears in brackets below it. 



The proposed solutions from Eq. [5] are substituted into 
Eq. [13] and solutions of the form 

5M= C k cos{{u + k)ujt) (14) 

fc=-oo 

are searched fo r, where v is a constant. CS determined from 
iFloauetl (|!883l Vs work that values outside the range ^ v ^ 
1/2 are unstable, and thus that v(xo) = and v(xo) — 1/2 
bracket the stable region. They found no solution in the 
case of v — 0, but solutions do exist for the case v = 1/2, 
discussed below. 

In order to get a finite number of non-trivial solutions, 
Eq. [TJ must be truncated after a finite number of terms. 
Following CS we take 

8M = Cfc cos (Jk + ujtj (15) 

which leaves us with six values of Ck to be determined. 

The main goal here is to solve for the six constants Ck- 
In order to do this, we substitute Eq. [15] and its derivatives 
into Eq. 1131 as well as the corresponding substitutions for N 
and D. From this point on, we diverge from the treatment of 
CS, as here we have used different methods to find this equa- 
tion's solutions. Here we used Maple 15 and Mathematica 8 
as tools for equation solving. 

1) Using Maple's COMBINE function, Eq. [13] was 
solved for one value of Ck- The solution revealed many co- 
sine terms with different frequencies, and some terms that 
were fully independent of the cosine. 



2) The three lowest frequencies of cosine (including the 
independent terms when present) were factored out of each 
individual Ck term. The terms relating to a single frequency 
were collected, yielding three separate equations. The cosine 
was then factored out, and the remainder of the equations 
set equal to zero. 

3) Steps 1 and 2 were repeated for each individual term 
of Ck ■ The result was 18 equations where there were 3 equa- 
tions for each Ck (each of the three equations representing a 
different frequency of cosine). Using these equations, a 6 x 6 
matrix results, where rows 1, 3 and 5 represent equations 
with k values of -3, -1, 1 and rows 2, 4, 6 represent equa- 
tions with k values of -2, 0, 2. 



" Su 





Sis 





S 1S 


■ 





S22 





8^4 





S26 


Ssi 





S33 





S35 








S42 





S44 





S46 


S51 





S53 





S55 








S S 2 





Se4 





Se6 . 



In other words, the columns are in increasing order from -3 
to 2, which demonstrates which equations contain which k 
values. 

We now solved the equations in sets of three for the 
values of Ck- In essence, we constructed equations from the 
matrix {e.g. Su + S13 + 1S15 = and so forth down the 
rows). Rows 1, 3, 5 were used to solve for C_3, C_i and 
Ci. Similarly, rows 2, 4, 6 were used to solved for C-2, Co 
and Ci- The two homogeneous sets of three equations were 
transformed to two non-homogeneous systems of order two 
with C-i and Co set equal to unity for mathematical con- 
venience. 

Once all six coefficients were determined, the values for 
x(t) and y(t) could be determined for specific sets of values 
of A, B and lj from Table [T] We then use the relations that 
u = - , and that r = ^/x 2 + y 2 to solve for u. 

Figure[2]shows two examples of a perturbed solution. In 
the left panel, a parametric plot in t of x versus y for parame- 
ters U — 2/3, C=0.1, and other values corresponding to line 
1 in Table [T] is shown. Here the initial value of y is taken to 
be 10 -4 to justify our assumption that it is small. The two 
solutions are so similar that the graphs overplot each other 
and cannot be distinguished. The right panel shows a much 
larger orbit based on the parameters in line 17 of Table [1] 
Here the Prendergast solution does not contain enough fre- 
quency information to completely reproduce the box orbit 
but captures some of the character of the true solution, such 
as the x and y amplitudes and period. 



4 THE APSIDAL ANGLE AND P-ELLIPSE 
ORBITS 

In this section, we turn our interest to the approximate so- 
lution of the radial orbital differential equation. We are in- 
terested primarily in precession of the apsidal angle, and so 
we will consider the problem now in terms of the anomaly 9 
rather than the time t. If one wished to determine the rela- 
tionship between these two variables, a Kepler-like equation 
would need to be solved. 
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CS considered the x and y equations of motion, but 
here we consider u = l/y/ x 2 + y 2 = 1/r with an eye to 
later using this result to determine the apsidal angle for the 
purely radial logarithmic potential. In this case, we examine 
the case where U ~ 1 and C -C 1. We start with the radial 
orbital differential equation, which is of the form, 



fu 
d9 2 



+ u = 



h 2 u 2 



f 



(17) 



where h is the angular momentum. We note that the force 
function /(-) is equal to —dV/dr, which can be obtained 
by differentiating Eq. [T] 

The Prendergast Method works very well for the so- 
lutions of the logarithmic potentials from Eq. [T] indicated 
in Section [2] for the x and y equations of motion and fur- 
ther elaborated in Section[3] However, the method does not 
seem to be well suited for purely radial logarithmic poten- 
tials with different initial conditions, and the solution does 
not agree with that obtained by pure numerical integration 
of the orbital differential equation. The p-ellipse a pproxi- 
mate solutions of the orbital equation, pioneered by IStruckl 
(|2006l ). is a much better way of not only deriving an accu- 
rate approximate solution to order e 2 (where e is the orbital 
eccentricity), but also obtaining the values for the apsidal 
precession to a high accuracy. We present a detailed analy- 
sis of the orbital equation for the radial logarithmic potential 
with or without the inclusion of the core scale length. In our 
analysis, the factor ^ C ^ 1 gives a measure of the core 
scale length (|Struckll2006T ). 

For the case of large orbits, or negligible core size soft- 
ening length, the equation of motion is given by 



uu + u = c 



(18) 



where the parameter c, in the notation of Struck, depends 
both on the constant scale mass M* and the core scale length 
s of the potential, the gravitational constant G, and the an- 
gular momentum h. The above equation is similar to EQ.fTTl 
and of the form 



u + u — 



Struck suggested an approximate solution of Eq. \T§\ 



u = - {1 + ecos [(1 - b) <j>]} 



1/2 



(19) 



(20) 



Here, p = a(l — e 2 ) is the semilatus rectum; a is the semi 
major axis (e. g. iMurrav fc Dermottl (jl999h : IValluri etajj 
|2005l )). e is the orbital eccentricity, and (1 — b) is the factor 
associated with the precession rate. Henceforth, for conve- 
nience, we set k = 1 — b in our analysis, and we will use 9 
instead of </> as was used by Struck. 

Struck finds that orbits obtained from a numerical inte- 
gration of the above differential equation look like precessing 
ellipses (p-ellipses) and considers the approximate solution 
given in Eq. [20] 

Substituting the solution of Eq. [20] into the differential 
equation, we find that 



2p 2 



2(1 + ecos kO) 



+ 



1 + e cos k6 



+ 



1 + e cos k6 
~2p2 



(21) 



k 2 1 - e 2 + 4(1 

UU + u = 



^)(1 + ecos kff) 



2p 2 



2(1 + ecos kO) 



— c 
(22) 



The LHS of Eq. [22] simplifies to 



— (1 + ecos k6) + i (1 + ecoske) + 
2p 2 p 2 



2p 2 



1 + e cos k8 



(23) 



Where ci = A- + A- = X and k 2 = (1 - bi) 2 = 2 are the 

2p- jj 1 

first approximations to k and c ([Struck 2006|). In a more ac- 
curate approximation to order e 2 , we find that the constant 
terms reduce to 



Ap 2 



1 



+ 



p- 



(24) 



Comparing next, the terms involving cos k8, we find 
that the coefficient of cos k9 is given by 

- (1 - e 2 ) - ft + ■ 
On simplification, one obtains 







(25) 



-=2| 1 + - + — 
1-he V 2 4 



(26) 



For e = 0, k 2 = 2 in accord with lStruckl (120061 ). 

In the case of non negligible core size, one has a similar 
though modified differential equation of the form 



(uu" + u 2 )(l + u 2 ) 



(27) 



The solution given in Eq. [20] upon substitution into the 
above differential equation leads to the expression 

k 2 , , „, e cos k6 k 2 
(1 + e cos k6) + H — + 



Ap 2 1 + e cos k6 



1 + 



1 + e cos kO 



(2 



which, upon comparison of terms independent of cos k8, sim 
plifies to 



1 \ k 2 



1 + 



1 



Ap 2 



e 



Comparing coefficients of cos k6, we get the more gen- 
eral dependence of k 2 . 



1 + 



1 + 4 



If terms of order e 2 are ignored, 



k z = 1 - bi = 2 



1 + 



1+4 



(30) 



(31) 



where bi is the first approximation of the precession factor. 
Furthermore, to order e 4 



k e 
Ap 2 



3 1 2 , 
-+—+—+ 

2 p 2 v 2 



fcV 
8p 2 



i + 4 



(32) 
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in agreement with Struck. 

It is interesting to note that the orbital differential equa- 
tion associated with apsidal precession is a special case of the 
generalized Burgers partial differential equations (GBE) and 
seems to characterize these equations similar to the way that 
Painleve e quations repres ent the Korteweg-de Vries type of 
equations (Sachdcv 1991). This variety of equations can be 
expressed as Eqs. 1331 and 1341 where f(x) and g(x) are suffi- 
ciently smooth arbitrary functions, a, e and c are real con- 
stants, and_^e_solutions of y are Euler-Painleve transcen- 
dents (|Kamkej|l943T ). 



yy" + ay' 2 + f(x)yy' + g{x)y 2 + ey + c = 



(33) 



In the case where f(x) and g(x) are constants, we have 
the Euler-Painleve equation 



yy" + ay' 2 + byy' + cy 2 + dy 1 a = 



(34) 



For a = 1, b = 0, the substitution y 
differential equation 



it 1 ' 2 leads to the 



lti" 1 
47 + f 



a u 

cu + d+ - = 

4 u 



(35) 



For a = 
results. 



1 the 4— terms cancel, and the following equation 



1 „ 
— u 

2 



cu + d = 



(36) 



It is important to observe that the orbital differential 
equation does not have the ay' 2 term in contrast to the 
GBE. This term contains terms of order e 2 and the correc- 
tion does not turn out to be significant. Hence, the p ellipse 
is a natural approximate solution of the generalized Burgers 
equations (GB E) and is an Euler - Painleve transcendant 
ijKamkdl 19431 ). 

As a rough estimate of the mean error in neglecting the 
^— term, we evaluate the following integrals that occur in 
the evaluation of this term. 



2fc 2 e 2 1 
16p e 2 



(1 + e cos k8) ? d9 - 



(1 + ecos key^dG 
(37) 



With the substitution kO = x we have 
(1 + ecos x) 5 dx 



h = 

and 
h = 



11 

7T k 



1 sin kn 
1+ 2 £ — 



8kn 



kn 
~2 



sin 2kn 



+ 



(38) 



1 1 

7T k 



dx 







1 sin kir 
1 e — ; 

2 kn 



(1 + e cos x) 2 
3 , 



16 



3 2 sm 2fc7T 

16 6 2/C7T + 



(39) 



Hence, we obtain for the difference of the two integrals, 



h-h 



sin kn 



e 

T 



2e sin2A;7r 



+ 



(40) 



kn 4 16 kn 
An approximate mean error (M.E.) due to the presence 



of the term ^— is 

Au 



M.E. = 



2fc 2 1 

167T p 



0.9e 
kn 



„2 „2 
e e 

T + 16 



(41) 



u 


e 




k 


e 


(rad) 


0.75 


0.1 


1 


81363 


i 


73222 


0.75 


0.5 


1 


87598 


l 


67464 


0.75 


0.9 


1 


99002 


l 


57867 


1 


0.1 


1 


73495 


i 


81077 


1 


0.5 


1 


81108 


i 


73465 


1 


0.9 


1 


98250 


l 


58466 


1.5 


0.1 


1 


62362 


l 


93493 


1.5 


0.5 


1 


78858 


i 


75647 


1.5 


0.9 


2 


17205 


i 


44637 


6 


0.1 


1 


43715 


2 


18599 


6 


0.5 


1 


53520 


2 


04637 


6 


0.9 


1 


98441 


1 


58314 



Table 2 . Some values of the apsidal angle from the p-ellipse pre- 
scription (C = 0) with varying eccentricity. 



Recalling p = a(l — e ) and taking k m 1.45 and e ~ 0.9 
2.9 1 



M.E. = 



16tt 0.19a 



We find that M.E. is ~ 1% for a = 1, e = 0.3 (Jfc = 1.79); 
M.E. increases with higer e and decreases with larger a. 
Interestingly, when kO — n or k9 = 0, 

12 

(43) 



Au 



showing that this correction term does not contribute to 
these angles, as shown below. 



(l-e)(l-e)2 (i_ e )i 



- + ^ ±£ T = -(l-e 2 ) + (l + e 2 ) 







(44) 



Next, we calculate the apsidal angle for the orbits in a 
logarithmic potential. The apsidal angle is the angle at the 
force centre between the smallest and largest apses, that is, 
between pericenter and apocenter. Hence, the behaviour of 
the logarithmic potential is similar to that of the n > 2 
power law potentials. Thus, there will always be a single 
minimum regardless of the value of the constant c. As c 
increases the location of the minimum simply shifts to larger 
x values. 

Only bound orbits are possible for this potential. As 
x — > oo, V(x) also approaches infinity due to the ln(x) term, 
so there is always an inner and an outer turning point no 
matter how large the total energy of the system. Stable cir- 
cular orbits are possible at the minimum of the effective 
potential. 

The approximate p-ellipse orbits are, on first appear- 
ance, only good to first order in e. However, Struck, in his 
thorough analysis, has shown that the orbital fits are excel- 
lent over several orbital periods. In fact, the value of k — 1 — 6 
which more accurately depends on e, is still fairly close to 
the more exact value; as partly due to the slow variation of 
c with e. The apsidal angle has been calculated for various 
values of e and is shown in Table f2] 

We now calculate the apsidal angle by using the Lam- 
bert W function, a function that is creating a renaissance in 
solving many interesting problems involving roots and limits 
of integration, as well as others. 

We begin by defining the energy E of an orbit through 
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the summation of its kinetic and potential energies: 

2" 



E = 



+ r 



+ V(r) 



(45) 



h is the angular momentum, and ^5 can be 
broken into 

Our main goal is to solve for 4| as this will provide us 
with an integrable function from which we can ultimately 
obtain a value for 6. 

Working in the regime where C C 1 and U ss 1, V(r) 
can be simplified further and we obtain the following 

d6_ _ h_ 

du 



\/2E + 2 In it 2 — v?A 



(46) 



where A — h 2 + 2C 2 and E = 0.5 + lnr c , where r c is the 
radius of the (c for 'circular') orbit. The value of u c = l/r c 
is taken here at values between 1 and 1.8, examining a range 
around the nominal value (E w 0, h = e~ 1,/2 ) of u c = e 1//2 « 
1.648. 

Where d6/du passes from positive to negative reveals 
the location of the apses, thus the limits of integration of 
Eq. [19] are its corresponding roots. We can solve for these 
roots by setting the denominator equal to zero, and manip- 
ulating it so that it can become solvable using the Lambert 
W function (I Valluri et alj|200ch . We start by reworking the 
denominator into the following form: 



1 2 U . 

mu A 



-E 



The roots are given by the expression 



-2WA- 



A 



(47) 



(48) 



where Wj represents the Lambert W function and j repre- 
sents the chosen branch. We solve for the two roots by using 
the -1 and the zero branches. 

Having the apocenter tm and pericenter r m distances in 
hand allows a determination of the orbit eccentricity through 



l + e 



(49) 



We note that solutions with imaginary eccentricity, which 
have two complex solutions which are conjugates of each 
other, would be manifested by a plunging of the orbi t into 
the force centre (|Hagiharalll93ll : IChandrasc khar 198^). 

Integrating Eq. [46] with the two roots as end points 
of the integral yields an answer that represents the apsidal 
angle for the particular orbit with a specific value of u c . 

Figure [3] shows the apsidal angle calculated by this 
method, for different values of C, E and h. 

Table [3] shows how values of u c ranging from 1 to 1.8 
yield si milar apsidal angles with va lues near For com- 
parison, iTouma fc Tremaind (|l9971 l use the epicyclic ap- 
proximation for near-circular orbits to determine that their 
g(a,y) (which is twice the apsidal angle as defined here) 
equals 2iv/y/2 = 4.44428 = 2 x 2.22144, a value close to 
the one arrived at here. However, a comparison with the 
numerically-derived result, also listed in Table [3] shows that 
the method proposed here is much more accurate: the two 
differ only in the fifth decimal place. As a comparison, we 
also show the apsidal angle for various values of small e using 



Lambert W 
approximation 



Numerical Difference 



1.0 


2.06310 


2.06300 


0.00010 


f.l 


2.07122 


2.07116 


0.00006 


1.2 


2.07868 


2.07862 


0.00006 


1.3 


2.08558 


2.08554 


0.00004 


1.4 


2.09201 


2.09200 


0.00001 


1.5 


2.09803 


2.09797 


0.00006 


1.6 


2.10368 


2.10361 


0.00007 


1.7 


2.10901 


2.10896 


0.00005 


1.8 


2.11405 


2.11397 


0.00008 



Table 3. The apsidal angle as calculated for different values of 
u c = l/rc- 



the p-ellipse approximation in the column labelled 'Numer- 
ical'. 

From Eq. 1261 we have 

-1/2 



fc = V2[l-| 



The apsidal angle is given by 

7T 7T / eN 1 / 2 7T 



2) 



V2 



32 



(50) 



(51) 



It is of interest to note that the roots can be found 
without any approximation for C in terms of the polylog 
function. For arbitrary C, one obtains from Eq. 0S]an exact 
expression 

InC 2 



1 



(52) 



for finding the roots of u . Now if we define k = In C — E 
and x = — c -l u -2 , Eq. [52] reduces to 

h 2 

kx + xln(l - x) = (53) 

* + ]*(!-*) = 4*x (M) 

k-Lhix) = ~ (55) 

Here Eq. [SS]is the fu nctional equa tion of the Polylogarithm 
Lii(x) = - ln(l - x) (|Lewinlll98ih and 



;■ 

C 2 



(56) 



5 GRAVITATIONAL LENSING 

The use of the Lambert W and the Polylogarithm functions 
to find the roots of equations such as Eq. 27] and [S3] may 
have wider applicability. For example, we can use a simi- 
lar approach to compute the deflection of a light ray by a 
logarith mic potential, u seful in the context of gravitationa l 



Cowlind Sill ; ISchutzlll99(]| : iBlundell et afll201Ch . 



lensing 

IZwickvl (1 19371 ) suggested that extragalactic nebulae offer a 
much better chance than stars for the observation of grav- 
itational lens effects. Zwicky's idea was that some of the 
massive and more concentrated nebulae may be expected 
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to deflect light by as much as half a minute of arc. Nebu- 
lae, in contrast to stars, possess apparent dimensions which 
are resolvable to very great distances . Zwicky was follow- 
ing up on the work of lEinsteml |l93rj ) on stars acting as a 
gravitational lens. According to Zwicky, observations on the 
deflection of light around nebulae may provide the most di- 
rect determination of nebular masses (|Smitblll936h . lZwickvl 
1 19371 ) estimated the probability of detecting nebular galax- 
ies which act as gravitational lenses and pointed out the 
possibility of ring shaped images, flux amplification and un- 
derstanding the large scale structure of the universe. The 
lensing equation can be generalized to three dimensions and 
cosmolog ical distances by corre ction of the redshift related 
distance ^Schneider et al.lll992l ). 

For arbitrary K one has the following expression to de- 
termine the roots in the case of light deflection for a loga- 
rithmic potential, 



V(r) = Kin 



(r 2 + C 2 ) 



(57) 



where K is a dimensionless constant and r = -. 

u 

For light deflection in the logarithmic potential consid- 
ered, the differential equation for the given logarithmic po- 
tential is of the form 



d 2 u 



+ u — K 



C l u 



-K 



dd 2 ' " "\1 + C 2 u 2 u) ' u(l + C 2 u 2 ) 
The DE for small values of C <C 1, reduces to 
d 2 u _ 2 -K 



d6 2 



KCu : 



(58) 



(59) 



In the relativistic formulation (|Hartle!l2003h the differ- 
ential equation for light deflection is 

2 1 

(60) 



_L 

i 2 + 2Mu 3 



u + u = 3M u 



(61) 



where b is the impact parameter. 

Assuming the photon is a non-relativistic particle that 
travels at speed c and it is fa r from all sources of gravita- 
tional attraction (|Hartlel [20031 ). we can determine the light 
deflection A9 produced by a logarithmic potential as 



A6 



du 



(62) 



/ l-tt 2 + A'ln(4 r +C* 2 ) 

Solving for the roots of the denominator, one obtains 

1 - u + K [In (1 + C V) - In u] = (63) 

Solving for K <C 1 by use of the Lambert W function, 
we get 



i K 



c 2 ) 



(64) 



where j denotes the branch of the Lambert W function. 

The defl ection angle is related to the Einstein angle 0e 
ljHartldl2003l ) which sets the characteristic angular scale for 
gravitational lensing phenomena. Gravitational lensing can 
be used to detect mass or energy in the universe, whether 
visible or not. Table|4]show the deflection angle for a range of 
values of K and C. Small K values produce small deflections, 
while smaller values of C produce larger ones, though with 



A 


C 


AO 


lO" 8 


1.0 

0.5 
0.0005 


0.00444 
0.00703 
0.00491 


0.005 


1.0 

0.5 
0.0005 


2256 33 
2861.51 
3202.15 


0.25 


1.0 

0.5 
0.0005 


69 637 4 
95 986.9 
109 909 


0.5 


1.0 
0.5 
0.0005 


104 102 
150 674 
173 731 


1 


1.0 
0.5 
0.0005 


142 042 
218 327 
251 802 


2 


1.0 

0.5 
0.0005 


178 640 
292 967 
334 185 



Table 4. The deflection angle A8 (in seconds of arc) as a function 
of k and C. 



a weaker dependence. Figure [4] shows the variation graphi- 
cally. 

Analogous calculations can be don e for time delay in 
light signals due to lensing galaxies l|Ohanian fc Ruffinil 
Il994f ). Intensity fluctuations caused by lumpy dark matter 
may provide direct observational existence for it. It is worth 
noting that the entire analysis can not only be done for the 
purely radial logarithmic potential, but also for the Eq.[T]of 
the logarithmic potential, by considering the x and y compo- 
nents separately as was done for the gravitational potential 
(IBourassa fc Kantowski|[l975l ). 



6 CONCLUSIONS 

We have revisited a nd e xpanded the work of 
IContopoulos fc Seimenisl (|l990l ) on orbits within a log- 
arithmic potential. We did a comprehensive review of the 
Prendergast method used by CS. We performed an analytic 
and numerical study of the matrix: C-3, C_2, C_i, Co, 
Ci, C2 that resulted in eighteen equations for the six 
coefficients for the orbital Fourier type series solution for 
values of U ranging from 0.1 to 1 that gave the unperturbed 
as well as perturbed solutions with better precision. The 
apsidal angle for the case of galactic orbits for a planar 
scale-free spherical logarithmic potential was obtained from 
the p-ellipse solution of the orbital differential equation 
and also the Lambert W . Both the Lambert W and the 
Polylogarithm functions may have applications in problems 
involving exponential and/or logarithmic potentials such as 
gravitational lensing. 

The Prendergast method, although not used as widely 
as others, has been quite useful in our analysis in Sections 
2 and 3., and is likely to prove useful in the study of many 
types of galactic potentials. 

Gravitational lensing can be used to detect mass in 
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the universe, whether dark or visible (|Hartlell2003l ; iNarlikarl 
2010). In general relativity, all energy curves spacetime, and 
a constant vacuum energy produces a detectable curvature. 
Gravity may prove a useful tool for detecting and studying 
dark energy. The lensing due to the gravitati onal field of a 
black hole of background stars and galaxies (|ThorneHl994T ) 
can be significant and the effects of a logarithmic potential 
warrant further study in this connection. 
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APPENDIX A: PERTURBED EQUATIONS 

The following are the equations for Equation ll6l 

Su = E e - U 2 E t (v~3) 2 lu 2 ; 
Sis = 0.5£ 7 - 0.5U 2 E 2 lu 2 (u - l) 2 - 0.5U 2 E 4 lu (y - 1) ; 
S 15 = 0.5E 8 - 0.5U 2 E 3 uj 2 (y + l) 2 - 0.5U 2 E 5 uj (v + 1); 

S22 = E 6 ~U 2 Et (v-2) 2 u) 2 ; 
S24 = 0.5E 7 - 0.5U 2 E 2 uj 2 u 2 - 0.5U 2 E 4 uv; 
S B6 = 0.5E 8 - 0.5U 2 E 3 w 2 (y + 2) 2 - 0.5U 2 E 5 uu (u + 2); 
S S i = 0.5E 7 - 0.5U 2 E 2 uj 2 (y - 3) 2 + 0.5U 2 E 4 uj (y - 3) ; 

Sss = E G - U 2 E t (u-l) 2 uj 2 ; 
S35 = 0.5E 7 - 0.5U 2 E 2 w 2 (y + l) 2 - 0.5U 2 E 4 w {v + 1); 
S42 = 0.5E 7 - 0.5U 2 E 2 w 2 (y - 2) 2 + 0.5£/ 2 £> (y-2); 

S44 = E a - U 2 E 1 u 2 v 2 ; 
S 4 o = 0.5E 7 - 0.5U 2 E 2 w 2 (y + 2) 2 - 0.5U 2 E 4 w {v + 2) ; 
Ssi = 0.5E 8 - 0.5U 2 E s u 2 (y - 3) 2 + 0.5U 2 E 5 uj (u-3); 
Sbs = 0.5^7 - 0.5U 2 E 2 w 2 (y - l) 2 + 0.5U 2 E 4 lu (u-1); 

S55 = E a - U 2 Et {v + 1) 2 oj 2 ; 
See = 0.5E 8 - 0.5U 2 E :i oj 2 (y - 2) 2 + 0.5U 2 E 5 lu {v - 2) ; 

Se 4 = 0.5E 7 - 0.5U 2 E 2 uj 2 v 2 + 0.W 2 E 4 uv; 

S 6e = E e ~U 2 E 1 {v + 2) 2 uj 2 . 

In S22 the bracketed term (y — 2) has been corrected 
from the original form with (y — 1) in CS. 



The corresponding E values are 

E, = C 2 (0.375B 4 + 3B 2 + 1) + 0.5A 2 (0.5B 2 + B + 1); 

E 2 = C 2 (3B 3 + 4B) + A 2 (0.375B 2 +B + 0.5) ; 

E s = C 2 (0.5B 4 + 3B 2 ) + 0.5A 2 (0.5B 2 + B) ; 

E 4 = 2cuB (A 2 + 2C 2 + C 2 B 2 ) + cjB 2 (0.5A 2 + 2C 2 B) - ljB 3 C 2 ; 

E 5 = loB 2 (A 2 + 2C 2 + C 2 B 2 ) + loB (A 2 + 4C 2 B) ; 

E 6 = U 2 uj 2 B (A 2 + 3A 2 B + 10C 2 B + 2.5C 2 B 3 ) + 0.75B 4 + 6B 2 + 2; 

E 7 = U 2 uj 2 B (2A 2 + 2.5A 2 B + 4C 2 + 13C* 2 B 2 ) + 6B 3 + 8B; 

E 8 = U 2 uj 2 B (A 2 - A 2 B + 2C 2 B + 2C 2 B 3 ) + B 4 + 6B 2 

In E4 the term with 0.5j4 2 has been corrected from its 
original form of 0.5^4 in CS. 
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Figure 1. A typical solution showing x versus time for the unperturbed solution based on the values of line 17 in Table[T] The solid line 
is the numerical solution from Maple 15, the dotted line is our approximation, and the dashed line is the difference between the two. 
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Figure 2. Two solutions showing x and y as a parametric curve in t. The left panel shows the Prendergast and real solutions based on 
the values of line 1 in Tabic [T] The two solutions plot on top of each other at this resolution and cannot be distinguished. The right 
panel shows a larger orbit, corresponding to the values of line 17 in Table [T] The solid curve is the numerical solution from Maple 15 
over a four cycles and the dotted line is the Prendergast approximation. 




Figure 3. The apsidal angle in radians as a function of the core radius C and the energy E. Four different angular momenta are presented 
in separate colours. Please see the text for more details. 
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Figure 4. The deflection angle in seconds of arc as a function of C and K. Please see the text for more details. 



